Localized Asymmetric Atomic Matter Waves in Two-Component Bose-Einstein 
Condensates Coupled with Two Photon Microwave Field 
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We investigate localized atomic matter waves in two-component Bose-Einstein condensates cou- 
pled by the two photon microwave field. Interestingly, the oscillations of localized atomic matter 
waves will gradually decay and finally become non-oscillating behavior even if existing coupling field. 
In particular, atom numbers occupied in two different hyperfine spin states will appear asymmetric 
occupations after some time evolution. 

PACS numbers: 03.75.-b,67.40.-w,39.25.+k 



The creation of dark solitons [l|, [l, 0] and bright soli- 
tons 0, 01 in trapped weakly interacting atomic gases 
has opened the road to understanding and controlling 
nonlinear properties of atomic matter waves. Although 
many nonlinear phenomena in Bose-Einstein condensates 
(BECs) may have their counterparts in nonlinear optics, 
unique manageability of the properties of BECs, new se- 
tups can be studied, which were not realized in nonlin- 
ear optics. This opens the promising perspective for nu- 
merous applications of the nonlinear matter- wave physics 
such as atom chip and quantum information processing 
on the nanometer scale. The s-wave scattering interac- 
tion in the BECs is the determining factor (attractive for 
bright solitons, repulsive for dark solitons) and attrac- 
tive interaction leads to condensates collapse in effective 
two- and three- dimensional cases if the atomic number 
exceeds a critical value, and to soliton formation in quasi 
one-dimensional (ID) traps. In the quasi- ID regime, two 
condensates in different hyperfine states reduce to cou- 
pled Gross-Pitaevskii equations (GPEs) [gI, In the 
system of optically coupled two-component BECs, some 
beautiful work [J S [13] has discussed Josephson-type 
density evolution process and pointed out the internal 
tunnelling effects of the two optically coupled conden- 
sates. 

In the present Letter, we consider two-component 
Bose-Einstein condensates coupled by two photon mi- 
crowave field. Using the soliton as staring spatial pat- 
tern formation, we study dynamic behavior of two com- 
ponents condensates under tuning atomic interactions. 
Surprisingly, the spatial shape of condensates presents 
very localized all the time even if changing the atomic 
interactions. During the transmission process, two local- 
ized atomic matter waves will collide each other many 
times and exchange energy greatly. In particular, oscil- 
lations of atom numbers of two components are very un- 
stable and will periodically decay to zero, and the speed 
of decay can also be controlled through tuning atomic 
interactions. Moreover, atom numbers occupied in two 
different hyperfine spin states will appear asymmetric oc- 
cupations in the end. 



Experimentally, the bright soliton have been success- 
fully produced from ^Li atoms in the internal atomic hy- 
perfine spin state \F = 1,^/ = 1) [4, 5]. Theoretically, 
we consider two different hyperfine spin state with at- 
tractive interactions denoted as and ^2, respectively. 
The ratio of the radial and the axial trapping frequency 
of harmonic potential obeys Ur/uJz ^ 1, for example, 
ujr/ujz = 5.1 X 10^) SO two hyperfine spin states 
are both in the transverse ground state. The conden- 
sate wave function is normalized as J (ir(|^ip + |^2p) = 
A^i + A^2 = where N is the total number of atom 
in two components. The attractive intra- and inter- 
atomic interactions are Uu = — 47rfi^|aii||^ip/m and 
Uij = —4:7TJi?\aij\\'^i\'^/m {i^j = 1,2), where — la^^l and 
— la^jl is the intra-species and inter- species s-wave scat- 
tering lengths, m is the atomic mass of the ^Li atom. 
In one-dimensional case the effective couplings Uu and 
Uij is represented by Uu = — 2^1^ \au\\'^i\'^ /mal and Uij 
= — 2fi^|aij||^ip/ma^ where = ^Jhjmijjj. [ll|. Af- 
ter tuning the strength of interaction to a sufficiently 
large negative value through use of atomic Feshbach res- 
onance [121, ll3| , the condensates can be set free along the 
waveguide, i.e. cj;^ = 0. The atomic transition between 
two-component BECs is induced by the two photon mi- 
crowave field with the effective Rabi frequency and 
a finite detune where the Vt and (5 are independent of 
both time and axial coordinate as experimental case [3] • 

In view of the difference of inter- and intra-atomic in- 
teractions, it is convenient to describe the two-component 
BECs coupled with the two photon microwave field in the 
following dimensionless coupled GPEs, 
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where z,j = 1,2, i ^ j and time t and coordinate 
z are measured in units 2/cj^ and = \/h/miJr, re- 
spectively, so, energy is in unit of huJr/2^ d = (|aii| — 
l^^il measure the difference of interaction strength, 
and / dz(\(j)i\^ + I02p) = 1 is normalization condition. 
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In order to compare with experimental case, some pa- 
rameters in equations ([T]) are chosen as cj^ = 27r x 625 
Hz (so, ar ~ 1.51 jim), \aiA = 0.16 nm (same as [4]). 



Since a^/|a^j| ~ 9437 
the number of ^ Li atom N 
4.N\aij\/ar ^ 2.5. 

In the special case of Ci = 0, i.e 



17, m 



so we choose 



([T]) have exactly two-soliton solutions 
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here i J = 1,2, i 7^ j and (9 = tan-i(f )/2, T : 
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Pj^ii)/{ki ~\~ k^){k^ ~i~ kj 
- 1^121^21)/ {ki + kl){k2 + /C2)l^i 
[(2iV|ai2|K)(a,a* + A/5*)][(^z 



A:*)].[(ai,/3i){A>:ii\A>:i2^} - {a2, p2){f^2i ^ ^^22}] 

In order to facilitate the understanding of the influ- 
ences of microwave field on solitons, it is convenient to 
obtain the asymptotic forms of solutions ([2]) of Eq. ([T]) 
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here superscripts in S'^^ {j^n = 1,2) denote soli- 
tons in limit t ±00, subscripts refer to the 
components and ^(j)'^'^) = {/cii^e*^^^sech(7^ii^ + 
i?i/2), k2Re'^-^sech{r]2R + R2/2)}, ((/>'-,(/>^+) 
= {k2Re''^'^sech[r]2R + (i?3 - i^i)/2], kiRc'^^^ sech[{r]iR ^ 
{Rs — i?i)/2]}, and the polarization vectors Aj 
and A]+are defined as {(yl^", A^"), (A^+, = 



[(27V|ai2|K)(|c^i 
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FIG. 1: Time development of the condensate densities (a) 
(b) |02|^ along the z axis for N = 6000, cjr = 27r x 
625 Hz, ttr ~ 1.51 /xm, 6 — 0.25u;r and Q — 0.75cc;r, 
ail — —0.24 nm, ai2 — —0.16 nm, 022 — —0.32 nm. The 
chosen parameters for initial solitons are ki — 0.313 + 0.2i, 
k2 0.313 - 0.2i, ai = f3i = (32 = 0.251 + 0.966i, 0^2 
0.966 + 0.25H. 
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|Ap)]M(c^i,A),(c^2,/32)} and 



{(A?-,Ar),(A{+,A^+)} = _c[(27V|ai2|K)(l«iP + 
IAP)^]{(«iM).(«2/a2)}{(Q^i^n - Q^2/^^i\/3i/^n - 

P2l^2l)^i^l'^12 - C^2l^22^Pl'^12 " P2l^22)} 

which c = (l/|A^i2p — l//^iiA^22)~2 , ai = 

{ki + /c*) [cos(7rz + 7r){ki - kj){a*aj + PtPj)]^- 

It is easy to see from the Eqs. (j3j) that the density of 
solitons are \S^^{z,t)\^ = cos^ ^+ l^^^^p sin^ ^ + 

(-l)^|A[^^||A^^^|sin2i9cos(2n + |^^^|^ n,lj = 

1,2, / 7^ j, where ip""^ = tan-^ (A^^^/A^^) - 
tan~^ (^2/^/^S)- ^^^^ special case, the reason 

of density oscillation is just because of the presence of 
time dependence part (-Ij^lAJ'^HA^^I sin2l9cos(2n + 



^^)0^^(z, t)^, where the frequency of oscillation is 2r 



and can be tuned easily through changing the effective 
Rabi frequency 1], a finite detune S and radial trapping 
frequency uOr- In particular, when 6 = n7r/2 {n e even) 
or r = 0, the oscillation will disappear. 
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FIG. 2: Time development of the condensate densities (a) 
(b) |02|^ along the z axis for N = 6000, cjr = 27r x 
625 Hz, ar ~ 1.51 /xm, 6 = 0.25u;r and Q = 0.75a;r, 
ail — —0.32 nm, ai2 — —0.16 nm, a22 — —0.24 nm. The 
chosen parameters for initial solitons are ki = 0.313 + 0.2z, 
k2 0.313 - 0.2i, ai = pi = p2 = 0.251 + 0.966i, a2 
0.966 + 0.25H. 



More interesting case is see the dynamics caused by 
changing the nonlinearities, i.e. atomic interactions. 
Experimentally, changing the nonlinearity can be real- 
ized by Feshbach resonance. In this system, changing 
the nonlinearity corresponds to change Ci in equations 
([1]). So, we will investigate pattern formation under 
changing Ci. We perform numerical simulations of Eq. 
([1]) based on discretizing it in both space and time us- 
ing the known boundary conditions, e.g., |(/>j (^min, ^)| = 



|(/>j(^max,^)| 



where 



-20 and z^ 



20. 
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Initial states we prepared are two soliton solutions ([2]) 
of Eq. ([1]). The mismatch of an with a^j in initial 
configuration will result in perturbations of the conden- 
sates. In our numerical simulations, the parameters 
are chosen as ki = 0.313 + 0.2z, /c2 = 0.313 — 0.2z, 
ai = pi = p2 = 0.251 + 0.966Z, a2 = 0.966 + 0.251z, 
S = 0.25Qr and Q = O.TbcOr- So, the frequency of oscilla- 
tions is 2r = 2[{6/ur)'^ + (n/cor)'^]^/'^ = 1.58 correspond- 
ing to 0.8 msec considering time units. 

Figure 1 shows image of initial conditions ([2]) prop- 
agating about 560 periods for |aii| = 1.5|ai2| = 0.24 
nm and \a22\ = 2|ai2| = 0.32 nm. It is interesting to 
point out that in our numerical simulations, spatial pat- 
tern formation all the time preserves soliton-like shaped 
i.e., localized atomic matter waves during transmission 
process. Although introduced important the mismatch 
of an with a^j, their influences on dynamics of initial 
solitons are only velocity and amplitude. That is to say, 
increasing the absolute value of atomic interactions will 
make the condensates more local instead of destroying 
localized density distributions. 

As is shown in Figurel, in the beginning, the densities 
of two localized atomic matter waves are oscillate peri- 
odically and later, oscillations will quickly damp out for 
each components as shown in Fig. 1(a) and Fig. 1(b), 
respectively. The total density of two components is well 
conserved during the transmission process, but the densi- 
ties of each components will change significantly. This is 
very incompatible with special case |aii| = |ai2| = 1^22 1- 

One of the most different phenomenon from special 
case is that two localized waves for each component no 
longer transmit separately and on the contrary, they 
transmit closer and closer. We can see that the two lo- 
calized waves collide for the first time at about 90 msec. 
After the collision, the relative density of one of waves 
is enhanced and the other is suppressed obviously. As 
shown in Figure 1, the density of component one corre- 
sponding to Fig. 1(a) is suppressed and the other corre- 
sponding to Fig. 1(b) is enhanced at the same time. So, 
the interaction between two waves can be used to trans- 
fer energy from one to the other. At about 190 msec, the 
two waves collide for the second time and oscillations are 
very weak now. Subsequently, the two waves will keep to 
collide faster and become closer and closer. 

It must be pointed out that whether the atomic inter- 
actions have the same values or not, the two localized 
atomic matter waves for component one are always com- 
pletely overlapped with the others for component two 
in configuration space. So, from the experimental point 
of view, observable phenomena are always the combina- 
tion effects of two components. The oscillations effect 
as shown in Fig. 1(a) and Fig. 1(b) can be observed 
through kicking BECs so that different hyperfine spin 
components will be separated in configuration space. 

In two different hyperfine spin components, two local- 
ized atomic matter waves are move closer and closer after 
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FIG. 3: Time development of the atom numbers (a) compo- 
nent 01, (b) component 02 for N = 6000, cjr = 27r x 625 Hz, 
ttr ~ 1.51 /xm, 6 — 0.25a;r and Q — 0.75a;r, an — —0.24 
nm, ai2 — —0.16 nm, a22 — —0.32 nm. The chosen parame- 
ters for initial solitons are ki = 0.313 + 0.2i, k2 = 0.313 — 0.2i, 
ai^ pi= p2 = 0.251 + 0.966i, a2 = 0.966 + 0.25H. 



a few collisions. So, now the questions are can we make 
two waves close faster and transfer energy more efficient? 
In Figure 1 , the atomic s-wave scattering lengths are cho- 
sen as an = —0.24 nm, ai2 = —0.16 nm, a22 = —0.32 
nm and the relation of absolute value of an and a22 is 
|<^ii| < |<^22|- So, it is naturally to investigate the time 
development of condensates in condition of inversing the 
relation of I an I and I a22 1, i.e. |aii| > |a22|- It is very sur- 
prising that two localized waves not only close faster but 
also transfer energy more efficient than former case. This 
is illustrated in Figure 2. The atomic s-wave scattering 
lengths are chosen as an = —0.32 nm, ai2 = —0.16 nm, 
a22 = —0.24 nm, and the other parameters are chose as 
the same as before. 

Fig. 2(a) and Fig. 2(b) show the time development 
of the condensate densities |0ip and |02p along the z 
axis, respectively. Due to changing the relation from 
|<^ii| < |<^22| to I an I > |a22|, the density of one of lo- 
calized waves has enhanced greatly and the other has 
suppressed extremely only after the second time collision. 
And the density of enhanced localized wave of each com- 
ponent will no longer increase or decrease during trans- 
mission process after the second time collision. Moveover, 
enhanced localized wave of each component will fixed in 
the position all the time and the weaken localized wave 
will move towards the boundary after the second time 
collision. The future applications of these effects can be 
existed in atom optics, atom transport and quantum in- 
formation. 

Another interesting problem is the effects of chang- 
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FIG. 4: Time development of the atom numbers (a) compo- 
nent 01, (b) component 02 for N = 6000, ojr = Stt x 625 Hz, 
ar ~ 1.51 /xm, 6 — 0.25a;r and Q — 0.75oOr, an — —0.32 
nm, ai2 — —0.16 nm, a22 — —0.24 nm. The chosen parame- 
ters for initial solitons are ki = 0.313 + 0.2i, /c2 = 0.313 — 0.2i, 
ai= pi= p2= 0.251 + 0.966i, a2 = 0.966 + 0.25H. 

ing atomic interactions on atom numbers of each compo- 
nents. From Eq. (|3]), it is easier to see that due to applied 
coupling field, atom numbers of each hyperfine compo- 
nent will exchange periodically. Naturally, the question 
is whether the persistent oscillations can still be observed 
or not when tuning intra- and inter- atomic interactions. 
Figure 3 and 4 show the time development of atom num- 
bers of each component for different atomic interactions. 
In Figure 3, as the system evolves, oscillations of each 
component both enhanced gradually in the beginning. 
At about 10 msec, the intensities of oscillations of each 
component both achieve maximum, and then begin to 
decrease, until at about 38 msec achieve minimum. At 
time interval of 38 msec to 80 msec^ the intensities of 
oscillations again change from maximum to minimum. 
Then, it is very surprising that after achieving maximum 
for the third time at about 115 msec, oscillations will 
gradually damp out and disappear finally. Moreover, it 
is interesting to notice that more atom number occupies 
the hyperfine spin state 02 as shown in Fig. 3(b) and less 
atom number occupies the other hyperfine spin state 0i 
as shown in Fig. 3(a) when oscillations disappear. The 
ratio of atom numbers of each component is 0.32 : 0.68. 

The same effects can be observed in Figure 4. When 
the relation of atomic interactions changes as shown in 
Figure 4, oscillations will experience more periods than 
former case before damp out. The ratio of atom numbers 



of each component is 0.35 : 0.65 and slightly different 
from former case. 

So, the conclusion is that atomic interactions have 
strong effects on exchange of atom numbers in this sys- 
tem. Through tuning atomic interactions using Feshbach 
resonance, we can control relative atom numbers occu- 
pied different hyperfine spin states and their oscillations. 
That is to say, we can control total magnetic moment 
of condensates using above mentioned effects. It is very 
fascinating to observe these effects experimentally. 

In this Letter, we have found a new type of wave - 
localized atomic matter waves in two-component BECs 
coupled with the two photon microwave field. By tuning 
the atomic interactions, localized atomic matter waves 
will experience a damped oscillations. During the trans- 
mission process, two localized atomic matter waves will 
collide many times and exchange energy significantly. In 
particular, the atom numbers occupied different hyper- 
fine spin states will appear asymmetric transition finally. 
Our result gives some fascinating effects on dynamics 
which are unobserved in real experiments so far. This 
work may be important for controllable nonlinear atom 
optics. 
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